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Abstract 

^ ' The spectral fluctuations of a quantum Hamiltonian system with time-reversal symmetry 

are studied in the semiclassical limit by using periodic-orbit theory. It is found that, if 
' long periodic orbits are hyperbolic and uniformly distributed in phase space, the spectral 

form factor K{t) agrees with the GOE prediction of random-matrix theory up to second 
Q ■ order included in the time r measured in units of the Heisenberg time (leading off-diagonal 

approximation). Our approach is based on the mechanism of periodic-orbit correlations 
discovered recently by Sieber and Richter By reformulating the theory of these authors 
Ph I in phase space, their result on the free motion on a Riemann surface with constant negative 

curvature is extended to general Hamiltonian hyperbolic systems with two degrees of freedom. 
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1 Introduction 



in 

• One of the fundamental characteristics of quantum systems with classical chaotic dynamics is 

■ the universahty of their spectral fluctuations. This universality and the agreement with the 

predictions of random-matrix theory (RMT) was first conjectured by Bohigas, Giannoni and 
Schmit (BGS) [2]. It has been later supported by numerical investigations on a great variety 
of systems [S]- However, the necessary and sufficient conditions on the underlying classical 
^ ', dynamics leading to such a universality in quantum spectral statistics are not known, and the 

origin of the success of RMT in clean chaotic systems is still subject to debate. 

In the semiclassical limit, where the BGS conjecture is expected to be valid, the Gutzwiller 
trace formula i4j expresses the density of states p{E) = X^„<5(£' — E^) of the quantum system 
I as a sum of a smooth part 'p{E) and an oscillating part. The latter is given by a sum posciE) = 

{TTh)~^ ^7 cos{S^/h — 7r^^/2) over all classical periodic orbits 7 of energy E {S^ and p-y 
are the action and the Maslov index of 7, and is an associated amplitude). The energy 
correlation function, 

and its Fourier transform K{t), the so-called form factor, are given by sums over pairs (7,7') 
of periodic orbits. Here r is the time measured in units of the Heisenberg time Tjj = 2iTh'p(E) 
{Th = 0{ti^~f) for systems with / degrees of freedom). The brackets denote an (e.g. Gaussian) 
energy average over an energy width W much larger than the mean level spacing /S.E = 'p{E)^^, 
but classically small, W <^ E, so that {p)e — 7>{E)- By neglecting the 'off-diagonal' terms, i.e., 
the contributions of pairs of distinct orbits modulo symmetries, Berry showed that the spectral 
fluctuations of classically chaotic systems agree in the limit /i — > with the RMT predictions 
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to first order in r (r <C 1). Two different approaclies have been proposed to support the BGS 
conjecture to ah orders in r in the semiclassical limit. The first one is based on a mapping 
between the parameter level dynamics and the dynamics of a gas of fictitious particles [3, SI- 
The second one uses field-theoretic and super symmetric methods and applies to systems with 
exponential decays of classical correlation functions 0. 

The link between spectral correlations and correlations among periodic orbits was first put 
forward in |S]. It was argued in this reference that the BGS conjecture implies some universality 
at the level of classical action correlations. Recently, Sieber and Richter ^ identified a general 
mechanism leading to correlations among periodic orbits in chaotic systems with two degrees 
of freedom having a time-reversal invariant dynamics. This has opened the route towards an 
understanding of the universality of spectral fluctuations based on periodic-orbit theory only. 
The crucial fact is that an orbit 7 having a self-intersection in configuration space with nearly 
antiparallel velocities is correlated with another orbit 7, having an avoided intersection instead 
of a self-intersection, which has almost the same action and amplitude. In two special systems, 
the free motion on a Riemann surface with constant negative curvature (Hadamard-Gutzwiller 
model) ^ and quantum graphs jS], the pairs (7,7) have been found to give a contribution 
1^2 (r) = — 2r^ to the semiclassical form factor. This result is in agreement with the Gaussian 
orthogonal ensemble (GOE) prediction of RMT, 



The first term Ki(t) = 2r is obtained by using Berry's diagonal approximation. 

The purpose of this work is to extend Sieber and Richter's result to general hyperbolic and 
ergodic two-dimensional Hamiltonian systems. Unlike in j^, our approach does not rely on the 
concepts of self- intersections and avoided intersections with nearly antiparallel velocities, but 
rather focus on what corresponds to such events in phase space, namely the existence of two 
stretches of the orbit (for both 7 and 7) which are almost time reverse of one another. It will 
be argued that working in phase space has a number of advantages and may allow for easier 
generalisations to periodically driven systems and to systems with / > 2 degrees of freedom. A 
similar approach is presented in [TO]; an alternative approach, based on a projection onto the 
configuration space as in |T, is presented in [TT] . 

In section |21 we state the main hypothesis on the classical dynamics used throughout this 
paper. After having briefiy recalled the main ingredients of the theory of Sieber and Richter in 
section |31 a characterisation of the orbit pairs (7,7) in the Poincare surface of section is given 
in section 0] The unstable and stable coordinates associated with a pair (7, 7) are introduced 
in the following section. The leading off-diagonal correction K2{t) = — 2r^ to the semiclassical 
form factor is derived in sectional Our conclusions are drawn in the last section. Some technical 
details are presented in two appendixes. 

2 Hyperbolic Hamiltonian systems 

We consider a particle moving in a Euclidean plane (/ = two degrees of freedom), with Hamil- 
tonian H{q,p) = H{q,—p) invariant under time-reversal symmetry. We assume the exis- 
tence of a compact two-dimensional Poincare surface of section S in the (four-dimensional) 
phase space F, contained in an energy shell H{q,p) = E and invariant under time reversal 
(TR) EI • Every classical orbit of energy E intersects S transversally. The classical dynamics 
can then be described by an area-preserving map (/) on S, together with a first-return time map 
X & Ti ^ tx ^ [0, cxd] (see J2]). In what follows, letters in normal and bold fonts are assigned 
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2r - rln(l + 2r) 



< r < 1 
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to the canonical coordinates x = {q,p) in E and to points x = {q,p) in F, respectively. It is 
convenient to use dimensionless q and p by measuring them in units of some reference length L 
and momentum P. The n-fold iterates of x by the map are denoted by x^ = 4>"' x, n € Z. They 
are the coordinates of the intersection points Xn of a phase-space trajectory with S, according 
to a given direction of traversal. The Euclidean distance between two points of coordinates x 
and y in S is denoted by |y — x|. If the system is a billiard {H{q,p) = /2M if q is inside 
a compact domain 17 C and +00 otherwise), S is the set of points {q,p) S T such that 
q is on the boundary dO, of the billiard, p is the momentum after the reflection on dO,, and 
\p\ = \/2ME. Then q is the arc length along dfl in units of the perimeter L, p is the momentum 
tangential to di} in units of V2ME, and tx is the length of the segment of straight line linking 
two consecutive reflection points, multiplied by the inverse velocity M/2E (see Fig. Due 

to the Hamiltonian nature of the dynamics, the linearised n-fold iterated map 
is symplectic. This means that it conserves the symplectic product 

Ax A Ax' = Ap Ag' - A(7 Ap' (3) 

for any two infinitesimal displacements Ax = (Aq, Ap) and Ax' = (Ag', Ap') in the tangent 
space T^S. 

The time reversal (TR) acts in the phase space F by changing the sign of the momentum, 
Tt '■ (q^p) ^ (Q) ~p)- Its action on x is given by an area-preserving self-inverse map T. When 
acting on an infinitesimal displacement Ax in T^E, the same symbol T refers to the linearised 
version of T (we avoid the cumbersome notation D^T, the meaning of T being clear from the 
context). In most cases, the exact map T is already linear and given by T : {q,p) {q, — p). 
The TR symmetry of the Hamiltonian implies (j)T = T(/>~^, i.e., {Tx)n = Tx-n- 

Some spatially symmetric systems in an external magnetic field have non-conventional TR 
symmetries, obtained by composing Tr with a canonical transformation associated with the spa- 
tial symmetry (SJ. The Sieber-Richter pairs (7, 7) of correlated orbits also exist in such systems, 
although they look different in configuration space JHI- By performing the canonical transfor- 
mation to redefine new coordinates {q,p) at the beginning, the TR becomes the conventional 
one. Therefore the analysis below also applies to systems with non-conventional TR symmetries. 

The normalised ^-invariant measure is the Liouville measure d/x(g,p) = dgdp/|S|, where |S| 
is the (dimensionless) area of S. Our main assumptions on the dynamical system (</>,S,/i) are 

(i) /u is ergodic; 

(ii) all Lyapunov exponents are different from zero on a set of points x of measure one (complete 
hyperbolicity); 

(iii) long periodic orbits are 'uniformly distributed in S'. 

Note that (i) and (ii) imply that the Lyapunov exponents itA^; are constant /i-almost every- 
where and equal to ±(A), with (A) > (the periodic points are notable exceptions of measure 
zero where this is wrong!). Examples of billiards satisfying (i-ii) are semi-dispersing billiards 
(if trajectories reflecting solely on the neutral part of dil. form a set of measure zero), the sta- 
dium and other Bunimovich billiards, the cardioid billiard, and the periodic Lorentz gas (see 
e.g. pi] and references therein). Assumption (iii), associated with ergodicity (i), means that an 
(appropriately weighted) average over periodic orbits with periods inside a given time window 
[T,T + 6T] can be replaced in the large-T limit by a phase-space average |15| I16j . Note that 
this statement, which is the precise content of (iii), does not concern individual periodic orbits 
but rather averages over many periodic orbits with large periods. We think that the statement 
can hold true even if some periodic orbits with arbitrary large periods are stable, if there are 
exponentially less such orbits than unstable orbits. 
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L = l 




Figure 1: Billiard qn and Pn = sin/?„ are the arc length along and the momentum 
tangential to d^l in dimensionless units. 

As is typically the case in billiards, the Poincare map (p or its derivatives may be singular 
on a closed set 5 C S of measure zero. For instance, if the boundary dQ is concave outward, 
(j) is discontinuous at a point xs = (qs^Ps) such that the trajectory between and the next 
reflection point is tangent to dfl at this point (see Fig.^. Let us denote by d{x, S) the Euclidean 
distance from x to S. We assume that 

(iv) S is 'not too big': fi{Bs^s) < Ci ^"^^ for any 5 > 0, with Bs^s = {x G S; |x— x^l < 5, G S} 
and o"! > 0; 

(v) the divergence of the derivatives of (p on S is at most algebraic, \d^<j)/dx"^ ...dx"''\ < 
C2d(x,5)-'^2(r-i)^ ^itj^ ^2 > 0. 

Here Ci > and C2 > 1 are constants of order one and the indices a,P = 1,2 refer to the q- and 
p-coordinates in E (x^ = q, x'^ = p). (iv-v) are standard mathematical assumptions on billiard 
maps [T7] . 



3 The theory of Sieber and Richter 

The starting point of Sieber and Richter is the semi-classical expression of the form factor, 

\ ^ \ T<{T-,+T^,)/2<T+ST I ^ 

The sum runs over all pairs of periodic orbits (7,7') such that the half-sum of their periods 
(T^ -|- T^i)/2 is in the time window [T, T -|- 5T] of width 5T <C Th- For isolated periodic orbits, 
= r^^ \ det(M-^^^ — where is the repetition (number of traversals) of 7 and M^^^ 
is the stability matrix of 7 for displacements perpendicular to the motion [3]. In order to work 
with a self-averaging form factor ^Hj, a time averaging over the window [T, T + 6T] (with, e.g., 
5T = h/W) has been performed in (jlj. Equivalently, K{t) can be defined as the truncated 
Fourier transform 



of the energy correlation function . Formula (HJ) gives the correct form factor for small enough 
times T only. It relates the quantum energy correlations to the classical action correlations Jj. 
Indeed, only orbits with correlated actions, differing by an amount of order h, can interfere 
constructively in Q. 

For fixed r = T/Th > 0, the sum Q) deals with orbits with very long periods as ^ — > 
(recall that Th = 0{h~^)). Such orbits have many self-intersections in q-space, some of them 
characterised by small angles e at the crossing point. As shown in the two loops at both 
sides of the crossing point can be slightly deformed in such a way that they form a neighbouring 
closed orbit in q-space, having an avoided crossing instead of a crossing (see Fig. 121). The two 
partner orbits 7 and 7 are almost time reverse of each other on one loop (right loop) and almost 
coincide on the other (left loop). Such a construction, which was supported in [T] by using the 
linearised dynamics, is in general possible in systems with TR symmetry and for small enough 
e only. Due to the hyperbolicity of 7, the two orbits come exponentially close to each other 
in q-space as one moves away from the crossing point q^. This means that the phase-space 
displacement perpendicular to the motion associated with 7 and 7 is almost (but not exactly) 
on the unstable manifold of 7 at Xi^c = (QcPcJ) whereas the displacement associated with 7 
and the TR of 7 is almost on the stable manifold of 7 at Xi^c [^^ ■ If a symbolic dynamics 
is available, the symbol sequence of 7 can be constructed from the symbol sequence of 7 in a 
simple way ^Hlj in the Markovian case, the TR symmetry implies that the partner sequence 
must not be pruned. Since the two orbits 7 and 7 have almost the same period and almost the 
same Lyapunov exponents, the amplitudes and are almost equal. Furthermore, it can be 
shown by using a winding number argument that = fi^ |lUllllj . In the Hadamard-Gutzwiller 
model, the difference 5S of the actions of 7 and 7 is given by 6S ~ E /\^^^ in the small e limit, 
where A^^^ is the positive Lyapunov exponent of the Hamiltonian flow pp. The main hypothesis 
of ^ is that, if the system has no other symmetries than TR, only the pairs (7,7) contribute to 
the leading off-diagonal correction K2{t) to the semiclassical form factor @ in the limit ^ ^ 0, 



The main task is to evaluate the right-hand side. This was performed up to now for the 
Hadamard-Gutzwiller model ^ and for quantum graphs The main difficulties arising in 
extending the theory of Sieber and Richter to other systems satisfying the hypothesis in the 
previous section are 

• the orbit 7 may have a family of correlated self-intersections, corresponding to one and 
the same partner orbit 7; this happens for instance in focusing billiards |111I2U| : care must 
be taken to avoid overcounting the pairs (7,7); 

• the specific property of the Hadamard-Gutzwiller model is that all orbits 7 have the same 
positive Lyapunov exponent A*^^-*; this clearly does not hold in generic systems; then the 
action difference 5S expressed in terms of e depends in general on 7; 

• the singularities X5 G 5 of the map (p affect the number of self-intersections with small 
crossing angles e and may even 'destroy' the partner orbit 7 if 7 approaches S too closely. 

We shall see in the following sections that working in the Poincare surface of section enables 
one to resolve all these difficulties. 




(6) 
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Figure 2: Pairs of correlated periodic orbits 7 (solid line) and 7 (dashed line) in configuration 
space for systems with conventional TR symmetry. The intersections of q-space with S are 
schematically represented by parallel vertical lines. 



4 The phase-space approach 

4.1 Orbits with two almost time-reverse parts 

As noted in 13 , if an orbit 7 has a self-intersection at with a small crossing angle e in 
configuration space, there are two phase-space points Xc^i = {q^,p^^) and x^j = {q^,Pcj) on 7 
which are nearly TR of one another, Xc^i ~ T^Xcj- Indeed, \p^^^ + Pcj\ — \Pc,i\\^\ is very small 
for |e| <C 1 (see Fig. |2]). There is in fact a part of 7 centred on Xc,i almost coinciding with the 
TR of another part of 7, centred on Xcj. The smaller the distance between T^Xcj and £Cc,i, 
the longer are these parts of orbit. There is therefore a family of points Xm of intersection of 7 
with the surface of section S, with coordinates Xm = 4>^x such that 

rx„_m PS , m = 0, ±1,±2, . . . (7) 

(see Fig. 121a)). The integer n is the time (for the map) separating the two centres x = xq and 
Xn of the two almost TR parts"*^ of 7. 

It turns out that the breaking of the linear approximation (LA) plays an essential role in 
the existence of a family {xm} with property 0. Indeed, we will show that, if the orbit 7 is 
unstable and n is large, the displacements 

— TXn—m (8) 

cannot be determined from Ax = Axq by using the LA for m > n/2. In order to make 
quantitative statements, and with the aim of transforming into a precise definition, we 
introduce a small real number Cx\ depending on x and on an integer t, the latter denoting the 
current time. Loosely speaking, ci*'' is the phase-space scale at which deviations from the LA 
after t iterations of a point y near x start becoming important. More precisely, this number is 
defined as the maximal distance \ym—Xm\ between the m-fold iterates of y and x, for an arbitrary 
y K, X and an arbitrary time m between and t, such that the final displacement yt — xt can be 
determined from the initial displacement y — xhy using the LA, yt — xt — Mx\y — x) (recall that 

^There is an analogy between the family of points {xm} and the family of vertices visited two times by an 
orbit in quantum graphs 0- 
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Mx is the linearised t-fold iterated map) ^. As the errors of the LA may accumulate at each 
iteration, the larger the time t, the smaller must be Cx^ . We shall see below that Cx^ decreases 
to zero like for large t if the map (j) is smooth. If cp is not smooth, a typical trajectory in E 
approaches a singularity point xs of <j) arbitrarily closely between times and t as t — > oo. As a 
result, Cx^ decreases to zero faster than at large t {cx^ even vanishes if x hits xs after m <t 
iterations, but such x form a set of measure zero). 

We can now define the time rriQ of breakdown of the LA for the displacements ((TJ as the 
largest integer such that 

lAx^l <4'"'') , m = 0,...,mo . (9) 

In other words, thq is equal to the largest integer m such that Axm — Ax. Similarly, going 
backward in time, we denote by rn^ the largest integer m such that Ax-m — Mx ™^Ax. In 
what follows, we say that the orbit 7 has two almost TR parts separated by n whenever @ holds 
true for a family {xm} of points of intersection of 7 with S, where Axm is defined by (jH)). The 
point xq is chosen among {xm} in such a way that |Axo| is minimum for m = 0. This condition 
fixes n. In order to simplify the notation, we shall drop the index for the coordinate xq of the 
centre point xq, writing x = xq and, similarly. Ax = Axq. Since we are interested in the limit 
|Axo| <C c^'^\ we always assume that mo and tuq are large (but much smaller than the period of 
7). If 7 is unstable, then |Axm| — |Mi"^^Ax| and |Ax_m| — \Mx Ax\ grow exponentially fast 
with m for large m with the same rate A-y, A-y = A^; > being the positive Lyapunov exponent of 
7 for the Poincare map. Moreover, the components of Ax in the stable and unstable directions 
are roughly the same, since, by assumption, |Axm| is minimum for m = 0. This implies that 
ttiq ^ mo- 

Let us first assume that n is large. The exponential growth of | Ax^l in the regime of validity 
m < rriQ of the LA has the following consequence. Let us look at the distance in configuration 
space in Fig. [21 between the point q^, moving on the lower branch of the right loop as one 
increases m (starting at m = 0), and its 'symmetric point' Qn-rm nioving backward in time on 
the upper branch of the same loop. After the time m = n/2, the two points on the lower and 
upper branches of the loop are exchanged. Thus, the distance between these two points cannot 
increase for m > n/2. In contrast, it must decrease exponentially as m approaches n, and come 
back to its initial value |q„ — QqI for m = n. It follows that the LA must break down before 
m = n/2, i.e., one has mo < n/2. A similar reasoning holds in phase space. We first note that 
the (n — m)-fold iterates of Txn and x are equal to Txm and Xn-m, respectively. Thanks to 
0, Ax„_m = —TAxm for any integer m. The equality |Ax„_m| = jTAxml would be violated 
if 2mo > n ^ 1, in view of the exponential growth of |Axm| predicted by the LA. As a result, 
for large n, the condition 

n > 2mo (10) 

must be fulfilled. 

Another situation arises when n is of order 1, n ^ mo- Then the above-mentioned arguments 
do not apply since, if m is of order 1, the unstable and stable components of Ax^ are of the same 
order and |Axm| does not necessarily increase with m. Since iQn-m ~ 9ml ^ \Txn-m — Xm\ *C 1 
for all times m between — mg and mo ^ n, the right loop in Fig. [21 consists of two almost parallel 
lines, connected by a small piece of line with length of order (/), (/) being the mean length of 
a trajectory between two consecutive intersections with S. This means that the orbit 7 has an 
almost self-retracing part in q-space, centred at Qq. 

more quantitative definition of ci'' and the precise meaning of 'yt —xt — M^\y — xy are given in section l^^ 

below. 
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Figure 3: (a) The two families {xm} and {xm} in the surface of section S (only few points are 
represented), (h) Magnification of (a) near x = xq. The N -periodic points x, x, Txn, and Txn 
pertain to 7 (filled circles), 7 (filled squares), the TR of ^ (empty circles), and the TR of j 
(empty squares). 

To conclude, we have shown that uq = 2mo has the meaning of a minimal time separating 
two distinct almost TR parts of 7 (i.e., excluding almost self-retracing parts). A similar result 
is obtained in jlUl lllj for continuous times. The continuous-time version of riQ is the minimal 
time To to close a loop in q-space introduced in pQ. In the present context, this time arises with 
the new interpretation of the breakdown of the LA. 

Let be the period of 7 for the Poincare map. If the family {x^} fulfils condition (jUJ, then 
the family of almost TR points {xn+m} also fulfils this condition, with n replaced by — n. 
This expresses the fact that, for a periodic orbit, the existence of a right loop in q-space implies 
the existence of a left loop (Fig. I^J. Setting y = Xn, one has Tyiy^n-m — Vm = —TAx^m and 

thus \TyN-n-m — Uml < ll^ll ci ^ for m = 0, . . . , . This indeed shows that if {x^} satisfies 
©, then this is also the case for {xn+m} with n replaced by — n, mo by tuq « m^, and niQ 
by m^' PS mo. 

The distinction made in the previous section between a self-intersecting orbit and an orbit 
with an avoided crossing in g-space is irrelevant in the surface of section S: both orbits have 
two parts which are almost TR of one another. In other words, they both have families of points 
{Xfn} and {xn+m} satisfying Q. Note that these two families can correspond in g-space with 
a family of self- intersections, as in the case of focusing billiards if self-intersections occur at 
conjugate points 12nj . or with one self-intersection (or one avoided crossing) only, as in the 
case of the Hadamard-Gutzwiller model pQ. 

4.2 The partner orbit 

We can now construct the partner orbit 7 described in section |31 in the surface of section S. Let 
7 be an unstable orbit of period A^ with two almost TR parts separated by n < A^. The orbit 7 
is defined by a A''-periodic point x = xq lying close to x = xq. This point is such that 



It can be checked in Fig. [3 that these properties indeed define the desired partner orbit. Note 
the symmetry of with respect to the exchange of x and x. By determining 6x = x — x as 




Txn-t — Xt\ for f = 0, . . 

xt — xt\ <^ 1 for t = n, . . 



n 



(11) 
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a power series in Ax, it is shown in appendix ^ that x has at most one partner point x. These 
arguments indicate moreover that x exists if |Ax| is 'sufficiently small' and the two almost TR 
parts of 7 are sufficiently far apart from a singularity point of (j). To first order in Ax, it is found 
in appendix IXI that 

Sx = i-x = T{l- Mi")Mi,^""^)"^(Mi") + T) Ax (12) 

in agreement with pp. The matrices Afi"^ and m!^^ appearing in ^Vl\ are the stability 
matrices of the right loop and of the TR of the left loop in Fig. |2I 

The partner point associated with Xm, — "t-q^ < m < mo, coincides with the m-fold iterate 
Xm = (f'^x of X. This can be seen by noting that Xm satisfies (|TT|) with x replaced by Xm and 
n replaced by n — 2m, as follows by combining (jllj) with @. Hence, by uniqueness, Xm is the 
partner point of Xm- It is not difficult to check this statement explicitly to lowest order in Ax on 
((T^ (see appendix IXjl . We conclude that the partner points of all points x^, — "Iq" < < tn^, 
belong to the same orbit 7. In other words, if |Ax| <C 1, there is a unique partner orbit 7 
associated with the whole family {xm}- If this family is almost self-retracing, i.e., if n <C mo, 
this orbit coincides with 7 itself, as already noted elsewhere 120]. Actually, then x = x satisfies 
(|1H) . hence 7 = 7 by uniqueness of the partner point (within the LA, this can be seen by replacing 
TAx = -Ax„ = -Mi"''Ax in ^] the identity Ax„ = Mi"^Ax follows from n < mo). By 
using a similar argument, one shows that the orbit 7' constructed from the family {Xn+m} is 
the TR of 7, as is immediately clear in Fig. |2 

4.3 A simple example: the baker's map 

The main advantage of the above-mentioned construction of the pairs (7,7) is that it works 
whatever the dimension of S (i.e., for systems with / > 2 degrees of freedom as well). Moreover, 
it applies to hyperbolic maps. It is instructive to exemplify this construction in the case of 
the baker's map. Then S is the unit square. It is convenient to equip S with the distance 
|x — x'l = max{|g — q'\, \p — p'\}. A point x = ((7,p) G S is in one-to-one correspondence with a 
bi-infinite sequence uj = ■ ■ ■ 00^2 ^-1 ■ ujoujiuj2 ■ ■ ■ , obtained from the binary decompositions of q 
and p {q = X];>o'^«-i and p = X^/>o'^-« 2~'), with binary symbols loi £ {0, 1}, I G Z. The 
map (j) acts on uj by shifting the point '.' one symbol to the right. The TR symmetry is the 
reflection with respect to the diagonal of the square, T : {q,p) ^ ip,Q)- This corresponds to 
reversing the order of the symbols of uj, i.e., T : uj ^ uj^ = ■ ■ ■ u)2 uji ujq . uj^i w_2 " " " • Periodic 
points are associated with sequences uj containing a finite word ujq - ■ ■ ujn-i, which repeats itself 
periodically; one usually writes the finite word only, keeping in mind that circular permutations 
of this word correspond to the same orbit. It is easy to see that the condition @ with ci™°^ = 2"'' 
is satisfied if a;„_/_i = uJi for any / = —s, . . . , mo -|- s — 1. Similarly, the condition | Ax_m| < 2~*, 
m = 0, . . . , rriQ , is satisfied if uJn-i-i = uji for any / = —rn^ — s, . . . , s — 1. This means that uj 
has the form 

X ^ uj = zIlZl.ZrRZI, (13) 

where Zl, Zr, L, and R are finite words containing {uiq + s), (mo -|- s), {N — 2mQ — 2s), and 
(n — 2mo — 2s) symbols, respectively. The symbol sequence of the partner point x is obtained 
by reversing time on R and leaving all other symbols unchanged, 

X < — > u = zIlZl. ZnR^Zji . (14) 

The inequality n > 2mo -|- 2s must be fulfilled in order that R is nonempty. In the opposite 
case, UJ has an almost self-retracing part ZlZjiZJ^ZJ^ and Co = uj. Similar pairs {uj,uj) of symbol 
sequences occur in the Hadamard-Gutzwiller model [T^ and in certain billiards ^201. The families 
{xm} and {xm} look like those in Fig. |2ta) after a rotation by an angle 7r/4. 
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5 Use of the unstable and stable coordinates 



To evaluate the leading off-diagonal correction K2{t) to the form factor, we shall first consider 
the second sum in @ over all partner orbits 7 of 7, for a fixed unstable periodic orbit 7, 
which will be assumed to be infinitely long and to cover densely and uniformly the surface of 
section. We will then argue in section El that one can replace the obtained result inside the sum 
over 7 in the limit T ^ 00. The sum over the partner orbits 7 of 7 is to be expressed as an 
integral over some continuous parameters characterising 7, chosen such that the action difference 
5S = S^y — S-y is a function of these parameters only. In configuration space, one may integrate 
over the crossing angle e P . It is argued in this section that a convenient choice of parameters in 
the surface of section is given by the unstable and stable coordinates of the small displacement 
Ax. The local coordinate system defined by the unstable and stable directions is singled out by 
the stretching and squeezing properties of the dynamics. These properties play a crucial role in 
the theory of Sieber and Richter, because they determine the time uiq of breakdown of the LA 
and the exponential smallness of the distances (fTT|) . 

5.1 The coordinate family Cx,r],£, 

Under the hyperbolicity assumption (ii), there are at almost all y € S two vectors eu{y) and 
es{y) tangent to the unstable and stable manifolds at y, which span the whole tangent space 
TyTi. These vectors can be found by means of the cocycle decomposition [T^ . 

Mf)e„(y)=AgeM(y^) , Mf )e,(2/) = A!™) e,(y„) , (15) 

where and A^^^ are the stretching and squeezing factors. Because My"^^ is symplectic, 

A^^^ = and the symplectic product e„(y) A es{y) is independent of y (see Jl]). The 

vectors eu,s{y) can be 'normalised' in such a way that this constant is equal to 1, 

eu{y)^es{y) = l. (16) 

The product of the norms of eu{ym) and es{ym) diverges as m — > ±00 if the angle between 
the unstable and stable directions at decreases to zero. Since the exponential growth of 
My^^Cuiy) at large m is (by definition) captured by the stretching factor, the divergence of 
|en,s(ym)| is smaller than exponential, ln\eu,siym,)\ = o(m) ^21- The notation f{m) = o{m), 
where / is an arbitrary function over integers, stands for f{m)/m as m ^ ±00. The 
stretching factor = A^ satisfy 

ln\A^J^^ = mXy + oim) (17) 

where Xy is the positive Lyapunov exponent at y. If y belongs to a periodic orbit 7 with period 
N, then eu,s{y) are the eigenvectors of the stability matrix My^^ of 7 and |Ay^^| = exp(A^A^). 
By invoking the TR symmetry, M^'^T = TMy ™\ Replacing this expression into ()15() . one 
finds eu,s{Ty) oc Te<j^u(y), with some i;^- invariant proportionality factors. By ergodicity, these 
factors are almost everywhere constant (y- independent). One can thus 'normalise' eu,s{Ty) in 
such a way that 

eu{Ty) = Tesiy) , es{Ty) = Teu{y) , A^^^ = A^ (ig) 

for almost all y G S. Note that this agrees with ((TB)) . since Tes{y) A Teu{y) = eu{y) A es{y). 
Two almost TR parts of an unstable orbit 7 can be parametrised by the family 

^x,'n,^ = {{rim, Cm) ; -"^o <m< mo} C (19) 
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Figure 4: Domain V^^a^u in the {ri,^)-plane (region marked by horizontal lines). The black points 
are the points {r]m,(,m) in ^x,r],^; ^ of them are contained in T^^y^ayU- 



of the unstable and stable coordinates {r]m,(,m) of the displacements Axm, 



Thanks to 



VmCm = r]C 



(20) 



(21) 



where -m^ <m< mo and fl = Vo-, = Co- The points {r]m,S,m) £ ^x,ri,e, are located on a 
hyperbole in the (r/, C)-plane (see Fig. 0]). 

In the case of the baker's map (section |4.3() . eu{y) and es{y) are independent of y and 
coincide with the unit vectors in the q- and p-directions. The stretching factors Ay"^ = 2"^ 
are also y-independent. The coordinates r]m and are the usual q- and p-coordinates of 



Cm 



Pn—m Qm 
Qn—m Pm 



2^{Pn-q) 
2-'^{qn-p). 



(22) 



5.2 Estimation of mo 

In the limit 77,^ — > 0, the time rriQ of breakdown of the LA depends logarithmically on the 
unstable coordinate ??, 

mo = -i^ + o(A-iln|7?|) . (23) 

/A'-y 

Indeed, thanks to hyperbolicity, jAx^l grows exponentially fast with m with the rate > 0, 
until it reaches, for m = rriQ, a value of the order of the phase-space scale ci"^°^ at which deviations 
from the LA start becoming important. Since |Axmol ~ l^mol' must have IndT/l/ci*""^) ~ 
-moXy. More precisely, we may approximate lAxyf^l by If/mll^M {xm)\ for m = ruQ and m = 
niQ + 1, making an exponentially small error for large mo (recall that |7y||e^(x)| ~ |C||es(x)|). By 
definition, |Axrrao| smaller than ci™°^ and |Axmp+i| is greater than c^°^^\ Then 1)23(1 follows 
from (fT7|). H21|). and ln|eu(xm)| = o(m). Note that the terms X^^\ Inci™"^! and A~^| lnc^x"°^^'^\ 
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have been neglected in H23|) . As stated in section H^Tl for a smooth map (j), ci"" decreases to zero 
hke mQ^^ as niQ oo, i.e., as — > 0. Therefore, | Inci*""^! is of order In mo = ©(Inj ln|r7||) and 
can be incorporated in the error term in H23|) . If (p has singularities on S, cl"") decreases to zero 
faster than mg ^ as r/ — > 0. In such case, it will be argued in section that formula (|23|) is not 
valid for all orbits 7. However, the right-hand side of H23|) always gives an upper bound on mo- 
Strictly speaking, the asymptotic behaviours H17|) and H23|) provide good approximations only if 
mo is close to a multiple of (or is much larger than) the period of 7. The physically relevant 
values of mo are, however, such that 1 <C mo <C A^. For such mo's, (|23|) should give nevertheless 
a reasonable approximation of the average value of mo (in fact it gives a good approximation 
of the inverse of the average of the inverse of mo). This average can be taken over all points x 
on 7 n S satisfying Q such that the unstable and stable coordinates of Ax = Tx„ — x are in 
small intervals [r/, r/ + dr/] and ^ + d.^], for an arbitrary integer n < A^/2 and some fixed r/, ^, 
di] <C <C 1, and d^ < |^| < 1. 
One shows similarly that 

ml = + o{X-Hn\^\) . (24) 

For the baker's map, in view of ()22() . 2™''|p„ — q\ < 2^* < 2"^'^^^\pn — q\, where we have chosen 
ci™"^ = 2~*. If mo ^ s, this yields mo — — In — In 2 and, similarly, mo" — — In Iq^^— p|/ In 2, 
in agreement with (|23|) and (|2H) . 



5.3 The probability of 'near-head-on return' 

To count the number of partner orbits of an orbit 7 with a very large period A^, one needs to 
know the probability to have two points on 7 which are nearly TR of one another. The aim of 
this subsection is to determine the (unnormalised) probability density P-y{r],^) associated with 
the unstable and stable coordinates of Axt = Txt+n — xt, for all pairs {xt,xt+n) of almost TR 
points on 7ns which do not pertain to an almost self-retracing family (i.e., such that n > 2mo). 
This density is defined through the number P-yirijS,) dr/d^ of points Xt on 7 n S such that the 
unstable and stable coordinates of Axt are in the infinitesimal intervals [ri,7] + dr]] and [C, C + d.^], 
for an arbitrary integer n between 2mQ{xt,rj) and A^/2. Let us recall that the partner orbits 7 
and 7' built from the two families {xm} and {xn+m}, separated from their almost TR families 
by the times n and N — n, respectively, are TR of one another (sectional). The two pairs 
(7,7) and (7,7') have thus identical contributions to the form factor © (the corresponding 
action differences 6S are clearly the same). This is why it suffices to consider only the orbits 7 
constructed from the family with the smaller time, n < A^/2. 
Let us define the infinitesimal parallelograms dS^, ,j ,f in S by 

dS^,^,5 = {y G S ; r/ < (y - x)n < r? + dr/ , e < (y - x), < e + de } (25) 
where (y — x)u,s are the unstable and stable coordinates of y — x. Then 

N-l N/2 

P^(r/,0d7?de= ^ J2 x(rx„+t e dS,,,^,^) (26) 

t=0 n=2'mo{xt,T]) 

where xi'P) equals 1 if the property V is true and otherwise. We shall assume here that the 
periodic orbit 7 covers densely and uniformly the surface of section S. If A^ ^ 4mo ^ 1, the 



12 



sum over n can then be replaced by a phase-space integral, giving 

P^^lr?, Od77de =J2i^Y~ J My)x{Ty G dS,„,,^) . (27) 

This integral is nothing but the area |Td5]a;j^^^^| = |dSa;j^^^^| of the parallelogram H25|) per unit 
area. By (fTO]) . it is equal to d7/d^/|Il|. In virtue of (^51) . 

F"%.0 = ^(iV + ii^) . (28) 

It is worth noting that the ergodic hypothesis implies the identity between ()26() and H27() for 
a set of points xq of measure one, and does not tell anything a priori about the points xq on 
periodic orbits, of measure zero in S. We shall argue below that, although (PB)) and ((77|) may 
differ for individual periodic orbits 7 which do not cover S uniformly, one can use the ergodic 
result ()28|) to calculate the form factor in the semiclassical limit. 

5.4 The domain V^a,!^ 

The density P'y{rj,(,) just defined overcounts the number of partner orbits 7 relevant for the form 
factor. Actually, a unique partner orbit 7 is associated with each family {xm] -ml<m< mo} 
(section 14. 2|) . whereas all points x^a belonging to the same family are counted separately in 
P-fiViO- To avoid overcounting, we define a domain T>^^a,u in the (??, ^)-plane M^, having the 
good property to contain, for any (77, ^) inside this domain, a fixed number v of elements {r]m^ £,m) 
in the family Cx^n,(- This integer u is independent of r/ and ^ (and thus of mo and ttiq) and 
is such that 1 <^ i' <^ N. Provided that this condition is fulfilled, the precise value of i' does 
not matter for the final result. Introducing also a small number a > controlling the maximal 
values of 1 77! and |^|, we define 

V^,a,u = {{^,0 e ; e-'^-' < 1§ < e''^-' ,\v^\<a^] . (29) 

is 



If (r?, ^) belongs to V^^a,!^, then \rj\ and |.^| are bounded by ae'^^''^'^. The domain 'D^^a,u 



represented in Fig. [l] For any (7?,^) G 1^j,a,i', it contains dx — v elements of the family Hx^-q^i- 
Actually, in view of (dH) and lf2T|). 

In ( 1^ ) = In ( 1^ ) - 2mA^ + o(m) , 1 <C |m| < min{mo, m^} . (30) 

Wlm\) \Vl\) 

By choosing a small enough, one has, thanks to mo + o(mo) > v for any (r/,^) G 'Dj^a,u 

(for instance, if ci™°^ = bmQ°' with 6, a > 0, one may choose a = be~'^'^^''^'^). The number of 
iVrm^m) in the family vC^.r;,^ which fulfil the first condition in H29|) is then equal to + o(z^). If 
(?7, ^) G 'D^^a,u, the second condition 1??™, < a is fulfilled, by (|2T|) . for all m between — mQ and 
mo, since it holds true for m = 0. Hence T>^^a,u H ^C^;^^^^ has d^,. = + o{v) elements. Note that, 
as already stressed in section 15.21 the use of the asymptotic behaviour H3(J() for 1 ^ mo ^ 
is in fact only justified if one is concerned with the average value of d^^, taken e.g. over all x 
on 7 satisfying © with unstable and stable coordinates of Ax in some intervals [r],r] + dry] and 
[.^, ^ + d^] for an arbitrary n < N/2. 

Let us define a new weighted probability density P^{r],^), in which the overcounting of 
partner orbits is compensated by a weight d"^ attributed to each event Txn+t £ dE^^^^^^g in (|26|) . 
By repeating the argument of the last subsection, one gets 

~.r., . ^ I fN , A 1 N f 4 1n|r7|\ 

^"'('^0 = E5;(y-2™o(-..''))p|-2ii>(" + ^) ■ '^^ 
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This density differs from (|2<S|) by a factor \/v. 
5.5 Action difference 

The main point in determining the action difference 5S = S^y — of the two orbits 7 and 7 
is to observe the following geometrical property of the partner points in the small |Ax| limit: 
X, X, Txn and Txn form a parallelogram, with sides parallel to eu,s{x) (see Fig. Efb)). It 
may be tempting to argue that, since, by (fTT|) . x must be exponentially close to the unstable 
manifold at x and the stable manifold at Tx„, this property follows straightforwardly from 
the continuity of the unstable and stable directions. However, some care must be taken here. 
Indeed, the unstable and stable directions vary notably inside the small region between the four 
A^-periodic points x, x, Txn, and Txn- This is due to the well-known intricate pattern built by 
the unstable and stable manifolds in the vicinity of heteroclinic points. We proceed as follows. 
Since Tx^ x = t] (x) +£^es{x), it suffices to show that, to lowest order in Ax, 

x-x = rieu{x) , Txn - X = £,es{x) . (32) 

The idea is to combine a stability analysis with the fact that eu{x) is nearly proportional to 
eu{x). For indeed, the orbits 7 and 7 look almost the same between times t = —{N — n) and 
t = 0. Therefore, their unstable directions must be almost parallel at x and x. Similarly, the 
TR of 7 is very close to 7 between t = and t = n, so that the stable directions at Tx„ and x 
must be almost parallel, es(Txn) oc es{x). 

To show (jSH), let us consider the unstable and stable coordinates {ip, C) of x — x, 

X - X = ip eu{x) + C es{x) . (33) 

In view of (|1H) . one may approximate Afi"^ by M^j?-^ and M^'^ by M^j^ if |Ax| <^ 1. By 
(|T^ . one has, to lowest order in Ax, 

-(1 - M^f )T{x - x) = (Mi") + T)Ax . (34) 

Here M^pp = Mj^-^Mj^ is the stability matrix of the TR of 7, with eigenvectors eu,s{Tx) 
and eigenvalues A^"^ such that |A^| = exp(A^A;y). By using (fTH|) . (pn|) . and and by neglecting 
terms smaller by a factor exp(— A^A^) or exp(— nA^) than the other terms, H34() can be rewritten 
as 

CA^e„(Tx)-Ve,(Tx) = Ai")r/e„(xn) + r/e,(rx) . (35) 

Hence, for n ^ 1 and (iV — n) ^ 1, C — and ipes{Tx) ~ —ij es{Tx). Replacing this result into 
and using ((TH|l , we arrive at the first equality in H32() . We now argue that the partner point 
of Tx is x„. This is already clear in Fig. |2 This can be shown by invoking the uniqueness of 
the partner point and by noting that the replacement of (x, x) by (Tx, x„) and of n by — n 
in leads to the exchange of the upper and lower lines, up to a TR. This replacement gives, 
by (jSJ, A(Tx) = TAx = ^ eu(Tx) + rj es{Tx). Then the second identity in H32() is a consequence 
of the first one (with the above-mentioned replacement), to which one applies the TR map T. 

The action difference 5S is determined to lowest order in Ax in appendix IbI It coincides 
with the symplectic area of the parallelogram (x, x, Tx„, Tx^), 

5S = {x — x) /\ {Txn — x) = 7/ ^ , (36) 

where we have chosen LP as the unit of action. It is clear that 6S is independent of the choice of 
the pair of partner points {xm, Xm), with — < m < niQ, as all these pairs {xm, Xm) correspond 
to the same orbit pair (7,7). Since r]m£,m is the only m- independent combination of rjm and S,rn 
of second order, the result (with an unknown prefactor) was thus to be expected. 
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6 Leading ofF-diagonal correction to the form factor 



6.1 The case of smooth maps 

The form factor © is, introducing a dimensionless Planck constant hcs = h/{LP), 

''■=M = ^J^( E *,d«A(,.?)exp(M)\ . (3T) 

The variables r] and ^ are integrated over the domain T>^^a,u defined in H29|) . As seen above, 
to avoid overcounting the partner orbits, one must use the weighted density P^{r],(^), related 
to the density -P7(?7, defined in section ESI by a factor l/v. Only partner orbits constructed 
from parts of 7 separated by n < N/2 from their almost TR parts are taken into account in 
these near-head-on-return densities, where N is the period of 7 for the map (p. The other partner 
orbits, corresponding to n > N/2, give the same contribution to the form factor (see section 1^^ . 
This contribution is taken into account by the factor 2 in 1)3 7() . 

The values of rj and ^ contributing significantly to the integral (|37|) are of order \/^eff • Thanks 
to (^5]) . no = 2mo is thus of the order of the Ehrenfest time A~^| ln?iefr|- For large periods, one 
has N ~ T/{ty) = r|S|/(27rncfr), where 

{ty) = I d^^{y) ty = mLP)-' I dy6{H{y) - E) = ^^"j^^^^^^ (38) 

is the mean first-return time. Therefore N ^ uq ^ 1 for the physically relevant values of t] in 
the semiclassical limit. This has also the important consequence that, for small but finite hcs, 
the values of the time T = t Th for which the theory of Sieber and Richter works are limited 
below by the Ehrenfest time 2Tq ~ 2{ty)no, since N must be bigger than 2no. 

We would now like to replace P'y{r],S^) by the ergodic result (pTT|) inside the sum (|37j) . To 
do this, one needs that long periodic orbits are uniformly distributed in phase space, in the 
sense explained in section El (see also |2II)- We shall assume here that this is the case, and that 
(|^T|) can indeed be used under the sum over periodic orbits (|n7|) in the limit T — > 00. A good 
indication supporting this assumption is given by Bowen's equidistribution theorem JH]: for 
any continuous function / on T, 

E e'"'-' '^-' T' dtf{x,{t))^ T.e-^'"''^-^ [ dfiEiy)f{y) (39) 

T<T-y<T+5T T<T-y<T+(5r 

( P) 

as T —)■ 00. The integral on the left-hand side is taken along 7, and X-y is the positive 
Lyapunov exponent of 7 for the Hamiltonian flow. The normalised microcanonical measure 
dfJ'Eiy) = S{H{y) — E) dy on the right-hand side is the product of the invariant measure fi 
and the Lebesgue measure along the orbit jl2j . 

l^dfiE{y)f{y) = l^dfi{y)F{y) , F{y) ^ j\t f{y{t)) . (40) 

Orbits 7 with multiple traversals > 2 have a negligible contribution in ()39() because they are 
exponentially less numerous than the orbits with = 1. One can thus replace exp(— A^^^ Ty) 
by the square amplitude A^^ in (jSH), 

E 4Ef(x„)~ NA'J dfi{y)F{y) , T ^ 00 . (41) 

T<T^<T+5T n=0 T<T-,<T+8T 



15 



To our knowledge, the sum rule (|H9|) has been proved rigorously for a restricted class of systems 
only, which includes uniformly hyperbolic systems fTB^ and the free motion on a Riemann surface 
with non-negative curvature 16 . Moreover, (|41|) cannot be applied directly to our problem, 
because x and mo in are discontinuous functions. We shall not pursue here in trying to 
motivate the above-mentioned assumption. Instead, we shall go ahead in determining K2{t). 
It would be interesting from a mathematical point of view to find general conditions on the 
dynamics implying our assumption. 

Replacing -By(?7, ^) by H31|) into the integral 

ha,u= [ d7?dCP,(77,e) expf^V (42) 



one obtains 



I-t,a,u = — < / — [ N + — — I sm 




A. 



7 / V 'ieflf 



sm 



2 



The first and third integrals can be computed with the help of the changes of variables ry' = 
rje^^-if^ and rj' = 7] e~'^^''^'^ , respectively. This yields 

The first term inside the brackets is equal to — vr -|- 0{hQg a~'^). The second one is a rapidly 
oscillating sine and gives rise to higher-order contributions in h^s after the energy average. 
Ignoring this oscillating term and the terms of order heg/a^, one gets 17,0,1/ = —T/Th- It should 
be stressed that this result is true only for very long periodic orbits which cover uniformly the 
whole surface of section S. It has been argued above that, although such a result is not true for 
all orbits 7, it can be used inside the sum over 7 in ((37j). This gives 

K2{t) = -^^ ^ A^^{l + 0{h,sa-')) . (45) 

H T<T-y<T+ST 

We can now invoke the Hannay-Ozorio de Almeida sum rule |21j . 

^ A2~r,r^oo, (46) 

T<T-y<T+5T 

to arrive at the announced result 

K2{t) = -2t^ , (47) 
valid in the limit h ^ 0, t = T/{27rhp{E)) fixed. 

6.2 The case of maps with singularities 

We have ignored so far the fact that (p or its derivatives may be singular on a closed set 5 C S of 
measure zero, as is typically the case in billiards ^7j- -^s stressed above, the term | Inci"*"^! 
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neglected in can be as large as mo if 7 approaches S too closely between times and mo- 
In such a case, it may a priori also happen that no partner orbit is associated with the family 
{xjn} (see appendix^. The aim of this section is to show that, under assumptions (iv) and (v) 
of section 121 the result ()47() is still valid. Indeed, we shall see that (|'23() and the action difference 
(|.S6|) are correct for all x outside a small subset of E. This subset turns out to be unimportant 
for K2{t) in view of its negligible measure. We will not discuss here the diffractive corrections 
to the semiclassical expression (jlj), which should a priori also be taken into account. 

Let us first estimate the phase-space scale Cx^ associated with the breakdown of the LA 
introduced in section [4.11 By invoking the cocycle property Mx^ = Mx]^-^ ■ ■ .Afij^AfiJ^ of the 
linearised map, it is easy to show that yt — xt is equal to 



M^\y-x) + \Y, f: MtlrUg^) [Mi")(y-x)]'[Mi™)(y-x)]%... , (48) 



where Mx^^ is the identity matrix. The displacement yt — xt can be determined from the initial 
displacement y — x by using the LA if the first term of the Taylor expansion (|48() is much greater 
than 
with 



than the subsequent (higher-order) terms. This is the case if \ym — Xm\ < Cx^ for < m < t 



(t) b ■ 
c\.' = - mm mm mm 

t m=0,...,t-l r>2 ai,...,ar=l,2 



. . . dx°'^ 



-l/(r-l) 

(49) 



A small fixed number 6 <C 1 controlling the error of the LA has been introduced. By assump- 
tion (v) of section 121 

4*^>-C2-^ min d{xm,Sy^ . (50) 



t m=0,...,t-l 



Let 6 > and 

mo-l 



r("io) 



U (l)-'^{Bs,s) , Bs,s = {x G S; |x - xs\ < 5 for some 3:5 G cSj . (51) 

m=0 

By (iv), it is possible to choose 5 such that the probability to find x in Bg^°\ 

mo-l 



m=0 

is very small. For instance, taking S = (b/mo) '^^ gives ^i{B^^s^) < < 1. Let us assume 

that X is not in B^"^°\ i.e., that the part of orbit between times and mo — 1 does not approach 
a singularity closer than by a distance S. Then, by (fSO)) . 

ct'^yC^H'^m,'^ , (53) 

with cr = 1 + o"2/o"i. Therefore, |lnci'"''^| is at most of order In mo = 0(ln | In |7/| |) as ?? ^ 
and can be incorporated in the error term in H23|). This reasoning shows that (\2'6\i can be used 
except if the centre point x of the family {xm} is in Bg^^^K 

If X is in B'f^°\ then mo may have a different behaviour for \ri\ <^ 1 than that given by 



H23|) . Anomalous behaviours due to singularities of the minimal time Tq to close a loop have 
been indeed observed in numerical simulations for the desymmetrized diamond billiard and the 
cardioid billiard in ^2 [201 • These numerical results show that non-periodic orbits satisfy (PS)) , 
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with replaced by the mean positive Lyapunov exponent (A), except those orbits approaching 
too closely a singularity. 

By using an expansion similar to (|48|) . one can show that the relative errors made by ap- 
proximating M^^ by M^)^ and M^^""^ by M^J""^ are small in the small \Ax\ limit if x is not 
in -Bj*^"^- The arguments of section 1^31 leading to the parallelogram {x,x,Txn,Txn) and to the 
action difference (|l-if)|) thus apply if x is not in i?^^). 

Let us now parallel the calculation of P-yirj,^) of sections and 15.41 Replacing the time 
average over t in (|^T|) by a phase-space average, 

P7S(7?,0 = N I dfi{x)^(^--2mo{x,7j)]y^ (54) 




-4 / _ ^d/z(x)mo(x,r?)+o(mo) 



with mo = — A^^ln|r/| and 6 = The integral in the second line gives a negligible 

contribution, as mo{x) < mo + o(mo) for any x (sectionEl and ^(5^^°^) < Cib. Thus 
is still given by with an error of order b. As 6 can be chosen arbitrarily small (in the limit 
^ ^ 0), it follows that K2{t) = — 2r^ for Poincare maps with singularities satisfying hypotheses 
(iv) and (v) of section [21 



7 Conclusion 

We have proposed a new method to calculate the contribution of the Sieber-Richter pairs of 
periodic orbits to the semiclassical form factor in chaotic systems with TR symmetry. Our basic 
assumption is the hyperbolicity of the classical dynamics. The method has been illustrated 
for Hamiltonian systems with two degrees of freedom. By assuming furthermore that long 
periodic orbits are uniformly distributed in phase space, the same leading off-diagonal correction 
K2{t) = — 2t^ as found in [T] for the Hadamard-Gutzwiller model has been obtained. This result 
is system independent and coincides with the GOE prediction to second order in the rescaled 
time r. One advantage of our method is its applicability to hyperbolic area-preserving maps, 
provided their invariant ergodic measure is the Lebesgue measure. This should allow one to treat 
the case of periodically driven systems. Moreover, the method is suitable to treat hyperbolic 
systems with more than two degrees of freedom /, for which the relevant periodic orbits do 
not in general have self-intersections in configuration space. A Sieber-Richter pair of orbits 
(7,7) is then parametrised by / — 1 unstable and / — 1 stable coordinates (r/^^\ . . . ,r]^^~^'^) 
and . . . ,^^-^~^^)- The time mg of breakdown of the linear approximation is given by the 
minimum of -ln|7?W|/A? over all i = 1, ...,/ — 1, where A^*^ is the ith positive Lyapunov 
exponent of 7. For Hamiltonian systems, the action difference 5S = — is given by the 
sum is 7- independent, whereas SS depends on the stability exponents of 7 in the 

approach of Sieber and Richter The evaluation of the integral (|37|) is more involved for / > 2 
than for / = 2 and will be the subject of future work. A second advantage of the phase-space 
approach is that it is canonically invariant and thus immediately applicable to systems with 
non-conventional time-reversal symmetries. A third advantage is, in our opinion, that orbits 
with crossings and avoided crossings in configuration space are treated here on equal footing. 

A further understanding of the universality of spectral fluctuations in classically chaotic 
systems may be gained by studying the contributions of the correlations between orbits with 
sei^era/ pairs of almost time-reverse parts ('multi-loop orbits') and their associated 'higher-order' 
partners. These contributions are expected to be of higher order in r. A first step in this direction 
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has been done recently for quantum graphs |22]- The phase-space approach presented in this 
work might be useful to tackle this problem. One would like to know if the RMT result @ can be 
reproduced in the semiclassical limit to all orders in r by looking at correlations between these 
partner orbits only, or if other types of correlations must be taken into account. An alternative 
way to study this problem is to investigate the impact of the partner orbits on the weighted 
action correlation function defined and studied in [51 123| . 

The periodic-orbit correlations discussed in this work have also remarkable consequences for 
transport in mesoscopic devices in the ballistic regime: they lead to weak-localization corrections 
to the conductance in agreement with RMT |24j . More generally, they should be of importance 
in any n-point correlation function of a clean chaotic system with time-reversal symmetry. 

Acknowledgements: I am grateful to S. Miiller for explaining his work on billiards, and to 
F. Haake and to one of the referees for their suggestions to improve the presentation of the 
results. I also thank P. Braun and S. Heusler for enlightening discussions and comments on this 
manuscript, J. Bellissard and M. Porter for their remarks on a first version of the paper, and 
H. Schulz-Baldes and W. Wang for interesting discussions. 

Note added. - When this work was mostly completed, I was informed that M. Turek and 
K. Richter were working on a similar approach. Drafts of the two papers were exchanged during 
the Minerva Meeting of Young Researchers in Dresden from January 29 to February 2, 2003. 



A Existence and uniqueness of a partner orbit 

We present in this appendix a general method, based on a Taylor expansion, to prove the 
existence and the uniqueness of the partner orbit. 

Let 7 be an orbit of period N with two almost TR parts separated by n < A^. Let x 
be the centre point of the family {xm] ~f^o < m < mo}. The partner orbit 7 is defined by 
an A^-periodic point x in the vicinity of x, called the partner point of x. This point fulfils 
property i.e., it is such that (i) \Txn-t — xt\ <C 1 between times t = and t = n, and 

(ii) \{Tx)t — {Tx)t\ <C 1 between t = and t = N — n. The small displacement 6x = x — x is 
obtained as a power series in Ax = Tx„ — x, 

00 

Sx" = (x - x)" = ^ [At'^] ^^...^^ Ax/^i • • • Ax^^ . (Al) 

r=l 

We use here the summation convention for the Greek indices /3i, . . . , /J^ = 1,2, referring to the 
q- and p-coordinates in S (x^ = q, x'^ = p). Let us stress that it is necessary to go beyond the 
linear approximation (term r = 1 in the series (jAll) ) to establish the existence of the partner 
orbit. Indeed, one must show that x is exactly A^-periodic, i.e., that xn = x to all orders in Ax. 

Let us assume that the map (p is smooth along 7 and its TR. We get the coefficients A^x^'^^ 
in HA1|) by expanding the final displacements as Taylor series in the initial ones for (i) the part 
of 7 between t = and t = n, and (ii) the part of the TR of 7 between t = and t = N — n. 
The identity xtv = x is then used to match the two results. More precisely, the computation is 
performed in four steps: (1) expand {Tx)]\f-n — (Tx)j\f-n in powers of Tx — Tx; (2) replace 5x 
by ()A1|) into this result; (3) expand Tx — Xn in powers of Txn — x and replace the series obtained 
in the previous step into this expansion and (4) identify each power of Ax. These manipulations 
lead for the linear order r = 1 to 

rA("'^) = M^") + T (A2) 



19 



with = 1 — Mx^^ ^\ Let us denote the partial derivatives {d^(j)^/dx^^ . . .dx^'')°' by 



[Mi*''')]^^ /3,' wi*^ t,r>l. For any 2 x 2 matrix C, we set [CA^x'''\...^p^ = Cl^[A^^''"\. 



The higher-order tensors A^''^\ r > 2, are obtained recursively through the formula 



s=2 riH |-»'3=»',ri>l 

^ /3r-rs + l-"/3r ' 



x[r4"-)]L_,...,^, (A3) 



with 



L a; Jai-.-a^ llsil ■ ■ ■ Sil J<5i---<5;L Ix iai---as^ 

1=1 siH \-si=s,Si>l 

We have assumed for simplicity that the TR map T on S is linear. 

It is worth noting that all tensors A^''^^ are obtained by inverting the same matrix If 
det ^i"^ ^ 0, then (jl^) reduces to (O and all A^'^'^'s are uniquely defined. Since tends 
to the stability matrix Mx^^ = Mx"'^ Mx^ of the unstable orbit 7 as I Ax I ^ 0, det# / Ofor 
sufficiently small |Ax|. This argument, however, does not suffice to show that is invertible 
for the physically relevant values of |Ax|, which are of order ^/TttH^ = Y^r|S|/A^ (section l6.ip . 
Another open mathematical problem concerns the convergence of the series ()A1|) . It can be 
expected that (|A1|) diverges when the orbit 7 approaches too closely a singularity xs & S 
between times —rriQ and ttt-q. Provided that is invertible and the series (|A1|) converges, the 
A'^-periodic point x exists and is unique. 

The above-mentioned construction is not restricted to the centre point x in the family 
{xm',—'mQ < m < mo}. Taking another point Xm in this family, one can as well construct 

its partner point (xm), by replacing x by Xm, n by (n — 2m), and Ax by Axm in ()A1|) . Let us 
show that, to linear order in Ax, (xm) is the m-fold iterate Xm of x. To lowest order in Ax, one 
finds 

£5 - X„ = ^ ikff )4«'1)Ax = x„ - x^ . (A5) 

The second equality is obtained by approximating 1 — D^"^ and M^™^ by M^^'* and 
respectively (see section 15. 5() . by using the cocycle property of the linearised maps, and by 
invoking the TR symmetry, which implies M^f^^T = T{My"^l^)~^ . It follows that (xm) = Xm 
belongs to the same partner orbit 7 as x. 

To conclude, we have given strong arguments in support of the existence of a unique partner 
orbit 7 associated with the family {xm',—fnQ < m < mo} if |Ax| <C 1 and the points in this 
family do not approach too closely a singularity of (p. 



B Action difference 

The action difference 5S = — between the two partner orbits 7 and 7 can be computed 
by considering separately the contributions 5Sji and 5Sl of the right loop (part of 7 between Qq 
and q„) and of the left loop (part between q„ and q^) in Fig. [3 dSn and SSl can be evaluated 
by means of the formula 

S{q„qf,E) - S{qi,qf, E) = {Pf + \ ^Pf) ■ ^Qf " (Pi + ^ ^Pi) ' ^<li + 0{\Sx{tt) , (Bl) 
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which gives the difference of action of two nearby trajectories q{t) and q{t) = q{t) + 5q{t) of 
energy E, initial positions q^ ^ q^ and final positions 7^ Qj. This formula is also valid for 
billiards. It is easily obtained by expanding the action difference up to second order in Sq^ and 
6qjr, and by using dS/dq^ = Pf and dS/dq^ = — Pj. In billiards, the momenta on the two 

trajectories have jumps a = p"^^ ~p\-cl ^ ~ P^rtl ~Ptc1 each reflection on the boundary 
do. (the sign — /+ refers to the values just prior/after the reflection). At first glance, a new term 
55'refl = — + cr) • <^Qrefi/2 should then be added to (|B1|) for each reflection point q^^^ on the 
unperturbed trajectory (such an additional term arises when writing the action difference as a 
sum of two contributions, corresponding to the two segments between q^ and q^^^ and between 
and qj). However, ^Srefi is of order |(^a3(t)p. Actually, 

5q,^^ = 5qT+-^KN + 0{5q^)=5qf--^KN + 0{5q^) , (B2) 

where 5q is the arc length on dO. between the two nearby reflection points q^^^y and Qi-efl, k is the 
curvature and T, N are the unit vectors tangent and normal to dVt at q^^^ (see Fig. The 
tangent vector T = T + 5qKN + 0{5q'^) at Qj-cfl appears in the last expression. Invoking the 
fact that cr and cr are perpendicular to the boundary, one gets (JSrefl = 0{5q^). 

Let us denote by x, x, x^ and Xn the points on the surface of section S with respective 
(g,p)-coordinates x,x,Xn and Xn- In the case of a billiard Q, these points are by definition 
associated with the values of the momenta just after a reflection on dQ. The corresponding 
points just before a reflection are denoted by the same letters with an added upper subscript 
(— ). The momentum jumps are denoted by <t = p — p^~\ with corresponding notation for p, 
p„ and p„. The action differences 5Sr and 6Sl are obtained by applying HB1|) with 

Xi = Tp Xyi , Xi = X , X J = Tp X , Xj = X12 

nf . — '1" o^. — nf nft J. — ^( ) nft j, — ™{ ) 

respectively. This yields 

2 = -(2p - p^) - p) • (g„ - g) - (-2p(,-) + p + p(^-)) • (q - gj (B3) 



2 = (2p(-^ + p(-^ - p^-)) • (q - q) - (2p„ + p„ - pj • (g„ - qj . (B4) 
A calculation without difficulties leads to 

26S = 6Sr + SSl = {x-x)a (Tr Xn - x) + (Tp x„ - Tr x^) A{x-Tr x^) 

-{a + (T)-{q-q)-&n- {qn " q) " • (q - q„) + ©(Ax^) . (B5) 

The F-symplectic product {y — x)f\{z — x) of two infinitesimal displacements {y — x) and {z — x) 
tangent to S at a;, with coordinates (y — x) and {z — x), reduces to the S-symplectic product 
{y — x) /\{z — x) given by Q (a choice of (g,j?)-coordinates in S with these properties is always 
possible, see Hence letters in bold font can be replaced by letters in normal font. The first 
term in the second line in ()B5|) is of third order in Ax by the above argument. One finds 

1 1 

dS = - (x - x) A [TXn -x) + - [TXn - TXn) A (x - TXn) 

-\ {{q - qn? -{q- qn?) f^nPn ' + 0{Ax^) , (B6) 

where k„ and Nn are the curvature and the normal vector of dil. at the point q^ of arc length 
Since x, x, Txn and Txn form a parallelogram to lowest order (see section x—Txn — x—Txn 
and the last term is of higher order in Ax. Therefore, (|B6() reduces to the canonical invariant 
expression Note that this result holds for any dimension of the phase space F. 
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